Strength of the Effective Coulomb Interaction at Metal and Insulator Surfaces 
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The effective on-site Coulomb interaction (Hubbard U) between localized electrons at crystal sur- 
faces is expected to be enhanced due to the reduced coordination number and reduced subsequent 
screening. By means of first principles calculations employing the constrained random-phase ap- 
proximation (cRPA) we show that this is indeed the case for simple metals and insulators but not 
necessarily for transition metals and insulators that exhibit pronounced surface states. In the latter 
case, the screening contribution from surface states as well as the influence of the band narrowing in- 
creases the electron polarization to such an extent as to overcompensate the decrease resulting from 
the reduced effective screening volume. The Hubbard U parameter is thus substantially reduced in 
some cases, e.g., by around 30% for the (100) surface of bcc Cr. 

PACS numbers: 71.15.-m, 71.28.-|-d, 71.10.Fd 



The effective on-site Coulomb interaction (Hubbard U) 
between localized electrons at surfaces of solids is ex- 
pected to be enhanced since the effective screening vol- 
ume of the surface is reduced with respect to the bulk. As 
a consequence, the electron polarization decreases at the 
surface, which reduces the effect of screening and gives 
rise to a larger U value. These arguments are underscored 
by interpolating between the Hubbard U values of an iso- 
lated atom and an atom in a bulk solid, the former being 
3 to 5 times larger than the latter jfj. Neither experi- 
mental nor theoretical works have been reported so far 
that would address the strength of the surface U param- 
eter explicitly. However, a large number of phenomena 
observed in solids indicates an enhancement of the U at 
surfaces. For instance, the metal-insulator transition at 
the surface of correlated materials [H , the appearance of 
magnetism at the surface of paramagnetic transition met- 
als [3| , and the enhanced exchange splitting at the surface 
of 3d ferromagnets 01 have been attributed to an increase 
of the correlation strength, which is defined by the ratio 
U /W , where W is the bandwidth. In the theoretical de- 
scription of surfaces, the U is usually assumed to be un- 
changed that the enhancement of correlation at 
the surface (S) with respect to the bulk (B) is provided by 
the effective band narrowing, i.e., W^g < W^. In princi- 
ple, depending on the relative values of the surface U and 
the bandwidth W , the correlation strength can further 
increase or decrease even below the bulk value. However, 
the latter case is considered to be unlikely because it is 
believed that U always increases at surfaces. In this Let- 
ter, we show by means of first-principles calculations that 
contrary to this conventional wisdom, this is not always 
the case. It decreases at many transition-metal (TM) 
surfaces and insulator surfaces with pronounced surface 
states, as a result of additional screening channels that 
open up due to surface-related changes in the electronic 
structure. 

Recently, the calculation of the Hubbard U parame- 



ter in solids from first principles has been addressed by 
several authors A number of different approaches 

have been proposed and applied to the bulk phase of var- 
ious classes of materials. However, the effective Coulomb 
interaction at a surface has been considered only within a 
model Hamiltonian framework [lij . so far. Reining and 



Del Sole [17| performed model calculations to account 
for the contribution of the surface states to the static 
electronic screening at the Si(lll) surface. The authors 
showed that the surface states give rise to a substantial 
enhancement of electron screening at the surface, reduc- 
ing the correlation strength. 

The aim of this Letter is to determine the strength 
of the on-site effective Coulomb interaction between lo- 
calized electrons at metal and insulator surfaces from 
first principles. To calculate the Hubbard U parame- 
ter we employ the constrained random-phase approxi- 
mation (cRPA) [ll| within the full-potential linearized 
augmentcd-plane-wave (FLAPW) method using maxi- 
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mally localized Wannicr functions (MLWFs) 
Our calculations show that the Hubbard U parameter 
is enhanced at simple metal (Na, Al) and most insu- 
lator (SrTiOs, NaCl) surfaces with respect to bulk as 
expected. However, the situation is different for TMs 
and insulators with pronounced surface states. For TM, 
both the interplay of the surface states and the effective 
band narrowing can give rise to a substantial reduction 
of the Hubbard J7, while for insulators only the surface 
states are responsible for the reduction of U . For the 
(100) surface of bcc Cr we obtain a 30% reduction of the 
static U . Moreover, the frequency dependence [/(w) is 
markedly different from that of the bulk. For bcc Cr we 
find that, starting from w = 0, the effective Coulomb in- 
teraction U{ijj) increases monotonically with frequency, 
and at about 2 eV it exceeds the bulk value, which re- 
mains basically constant in this interval. 

We model the metal and insulator surfaces with slabs 
of 11 atomic layers. Such slabs form a superlattice with 
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20 A of vacuum separating them, with each slab pos- 
sessing two (100) or (110) symmetric surfaces. To dis- 
cuss the chemical trends, additional modifications such 
as surface reconstruction and surface relaxation are at 
first not taken into account. We consider the 3d TMs 
in their respective ground-state crystal structures except 
Sc, Ti, and Co (Mn), which are treated in the fee (bcc) 
structure. The ground-state calculations are carried out 
using the FLAPW method as implemented in the FLEUR 
code [3l with the GGA exchange-correlation potential 
as parameterized by Perdew et al. [2^ . The MLWFs are 
constructed with the Wannier90 code [21, 22|. The effec- 
tive Coulomb potential is calculated within the recently 
developed cRPA method ll| implemented in the SPEX 
code [i^ (for further technical details see Refs.[l3 and 
H)- 



(a) Simple metal [fee AI (3p)] 



The cRPA approach offers an efficient way to calcu- 
late the effective Coulomb interaction U and allows to 
determine individual Coulomb matrix elements, e.g., on- 
site, off-site, intra-orbital, inter-orbital, and exchange as 
well as their frequency dependence. The basic idea be- 
hind the cRPA [llj is to define an effective interaction U 
between the localized electrons by restricting the screen- 
ing processes to those that arc not explicitly treated in 
the effective model Hamiltonian. To this end, the full 
RPA polarization matrix P is divided into P = P\ + P^, 
where Pi includes only transitions between the localized 
states, for which the Hubbard U is to be calculated, and 
Pj. is the remainder. Thus, the localized states are largely 
eliminated, and the screening is dominated by itinerant 
s and p states, which are well described by LDA and 
GGA. The cRPA is, therefore, - even with these stan- 
darad xc potentials - considered a reliable approach to 
calculate the Hubbard U parameter and its frequency 
dependence 2^, 2^. Then, the frequency-dependent ef- 
fective Coulomb interaction is given schematically by the 
matrix equation U{uj) = [1 — vP^{u!)]~^v, where v is the 
bare Coulomb interaction. The static limit of the average 
diagonal matrix element oiU{oj 0) represented in a lo- 
cal basis can then be regarded as the Hubbard U parame- 
ter. The matrix elements of the effective Coulomb poten- 
tial U in the MLWF basis are given by C/Rnin3;n4„2 ('^) = 
IKM^)wnsRir)Uir,r';uj)w*^^i^{r')wnMr^') d'^r <fir' , 
where w„R(r) is the MLWF at site R with orbital index 
n and C/(r,r';a;) is calculated within the cRPA. Wc de- 
fine the average on-site diagonal (direct intra-orbital) and 
off-diagonal (exchange inter-orbital) matrix elements of 
the effective Coulomb potential asU = L^^ Yin U-R.nn:nn 
and J = [L{L - 1)]"^ I]m,n(m5^«) ^R™";"™' where L is 
the number of localized orbitals, i.e., one, three, and five 
for ,s, p, and d states, respectively. 

In Fig. [T] we present the calculated layer dependence of 
the Hubbard U for three different systems: (a) for the Zp 
states of the simple metal fee Al, (b) the 2>d states of the 
transition metal bcc Cr, and (c) the Cl-3p states of the 
insulator NaCl in the rock-salt (RS) structure. As seen in 
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FIG. 1: (Color online) (a) Left panel: layer dependence of the 
Hubbard U parameter for fee Al. The S denotes the surface 
layer in the slab model. Right panel: Total density of states 
(DOS) for the (100) and (110) surface of fee Al. For compar- 
ison, the bulk DOS (shaded area) is included; (b) and (c) the 
same for bcc Cr and rock-salt NaCl. 



TABLE I: Bulk and surface Hubbard U and J values for 
simple metals and insulators. The corresponding orbitals for 
which the U and J are calculated are given in parenthesis. 
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ulk 


100 Surface 


110 Surface 




Ub 


Jb 


Us 


Js 
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Js 


Na [3s] 


1.39 




1.50 




1.47 




Al [Sp] 


2.18 


0.36 


2.46 


0.39 


2.36 


0.38 


MgO [0-2p] 


7.10 


0.63 


7.23 


0.62 


6.38 


0.59 


NaCl [C\-3p] 


7.13 


0.55 


7.60 


0.54 


7.53 


0.55 


SrTiOs [Ti-3rf] 


3.34 


0.37 


3.62 


0.39 






SrTiOs [0-2p] 


4.42 


0.56 


4.79 


0.55 







cases (a) and (c), the Hubbard U increases from the mid- 
dle layer, where it is close to the bulk value, to the surface 
layer as expected. However, we find an unexpected be- 
havior in the case of bcc Cr, where U is substantially 
reduced at the surface. This reduction is 30% for the 
open (100) and 20% for the (110) surface. For simplicity 
- the magnetism of bcc Cr is quite complicated - we only 
consider the non-magnetic state here. Furthermore, the 
layer dependence of U is quite different in metals and in- 
sulators. Because of the short screening length in metals, 
the Hubbard U quickly assumes the bulk value in the for- 
mer, as we go from the surface toward the middle of the 
slab, while in the latter the laycr-by-layer convergence 
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to the bulk value is much slower. On the other hand, 
the surface J values, which are listed in TablcUfor vari- 
ous materials, only differ slightly from the corresponding 
bulk values. 

We find that the surface electronic polarization and, as 
a consequence, the Hubbard U parameter is determined 
by two competing effects: (i) the so-called dimensionality 
effect, which is due to the reduced coordination number 
and hence the decrease of the effective screening volume 
at the surface region. From the point of view of clas- 
sical electrostatics, this effect reduces the electronic po- 
larization at the surface leading to larger U values, (ii) 
Electronic structure effects, i.e., the appearance of sur- 
face states and the effective band narrowing. This second 
effect gives rise to an enhancement of the electronic po- 
larization and, hence, to a decrease of U . Depending on 
the strength of the two competing effects, the effective 
Coulomb interaction at the surface can be enhanced as 
well as reduced with respect to the bulk value. Quali- 
tative information on the influence of the surface elec- 
tronic structure on the Hubbard U parameter, leading to 
the second effect, can be deduced from the polarization 
function, which is given by 
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P(r,r';w) = 2 ^ ^ ^„(r)C.' (r)C.(r')^m' (r') 



where '!/'in(r) and £,„ are Kohn-Sham eigenfunctions and 
eigenvalues, respectively, 5 is a positive infinitesimal, and 
Amm' = Em' " Em- For zero frequency (w = 0) the main 
contribution to the polarization function comes from the 
states around the Fermi energy. The smaller the energy 
difference between occupied and unoccupied states the 
larger the contribution. The effective band narrowing in 
TMs tends to reduce A^,^, with respect to Aj^,^, , which 
has the effect of increasing the polarization. Additionally, 
the presence of surface states close to the Fermi level at 
the TM bcc (100) surfaces make Af,^^, effectively smaller 
resulting in a more efficient electronic polarization and, 
as a result, in substantially reduced U and U values as 
shown in Fig. [5] for the id series, where U stands for 
the fully screened Coulomb interaction within the RPA. 
As seen in Fig.[Tl the density of states (DOS) of the 
(100) surface of bcc Cr exhibits a pronounced peak that 
is caused by a surface state just below the Fermi level, 
which is also found in previous first-principles calcula- 
tions [23| • This surface state contributes substantially to 
the polarization function and is mainly responsible for 
the 10% stronger reduction of the U value compared to 
the corresponding value for the (110) surface, where such 
a peak in the DOS is missing. On the other hand, in sim- 
ple metals the surface electronic structure turns out to 
be very similar to that of the bulk so that the dimension- 
ality effect (i) wins over the electron-structure effect (ii) , 
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FIG. 2: (Color online) (a) Partially screened {Ub) and fully 
screened (Ub) Coulomb interaction for the bulk 3d TM series 
in the non-magnetic state. With open diamonds we show 
the U for ferromagnetic Fe, Co, and Ni. (b) The difference 
between the surface and bulk U values. Open [filled] diamonds 
show Us — Ub for the (100) [(110)] surface of ferromagnetic 
Fe, Co, and Ni. (c) Bare Coulomb interaction V for bulk and 
surfaces, (d) The same as (b) for the fully screened Coulomb 
interaction U. 



giving rise to enhanced Hubbard U parameters. 

For most other 3d metals, the Hubbard U parameter 
is reduced at the surface, too, but less than in the case 
of bcc Cr, and for the late TMs the surface U exceeds 
that of the bulk U. In Fig.[5Ja) we present the partially 
screened Ub and fully screened Ub values for the bulk 
phases of the 3d TMs in the non-magnetic state. Results 
for the ferromagnetic ground states of Fe, Co, and Ni 
are also included for comparison. The difference between 
the surface and bulk U and U values, i.e.. Us — Ub and 
Us — Ub, is presented in Fig.[2jb) and (d) for the (100) 
and (110) surface. As seen, from Sc to Fe the Hubbard 
Us is reduced at both surfaces {Us — Ub < 0), and Co 
is at the border, in which the Hubbard U assumes sim- 
ilar values in the bulk and at the surface. Only at the 
Ni and Cu surface the U is slightly larger. As for the 
difference of the fully screened Coulomb interaction, i.e.. 
Us ~ Ub, we obtain a qualitatively similar behavior, but 
the relative reduction of U with respect to the bulk value 
is significantly larger. This is attributed to the fact that, 
in contrast to the Hubbard U, screening effects stemming 
from 3(i — > 3d transitions contribute to the effective in- 
teraction, too. At the surface these virtual transitions 
take place within the surface states, which leads to very 
effective screening effects that give rise to the observed 
reduction of the U. It is important to point out that 
the variations of U seen in Fig. [2] are not caused by dif- 
ferent spreads of the Wannier function across the series, 
as can be seen from the bare Coulomb interaction V, 
which is presented in Fig.[2][c). As seen, bulk and surface 
V values monotonically increase from Sc to Cu, with no 
apparent difference between bulk and surface values. We 
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FIG. 3: (Color online) (a) Frequency dependence of the bulk 
and surface Hubbard U parameter for bcc Cr. In the inset we 
expand the low frequency region, (b) The same for fee Ni. 



note that the Hubbard U values depend only little on 
the bare Coulomb interaction V because in metals we 
are in the strong coupling limit, i.e., f|-Pi| ^ 1, and thus 
U ~ — As for the surface Us of the ferromagnets Fe, 
Co, and Ni, the same discussion holds. So far, we have 
not taken into account the surface relaxation, which is 
usually small (< 5%) for most of the 3d TMs [H, 
Only the bcc V (100) surface possesses a sizeable inward 
relaxation (~ 11%), which results in a small change (~ 
5%) of the calculated surface U and U values. 

In contrast to metals, the Hubbard U at insulator sur- 
faces is much more strongly affected by the presence 
of surface states than the band narrowing. For exam- 
ple, the NaCl (100) surface docs not exhibit any surface 
states , and the U is enhanced. The slight changes in 
the electronic structure of the (110) surface (the gap is 
smaller, see Fig.[T]) reduces the surface U only very little 
so that it is still larger than the bulk value. The situation 
is similar for the SrTiOa surface, for which the results are 
presented in TableHl In MgO, on the other hand, both 
(100) and (110) surfaces exhibit surface states, which 
effectively reduce the band gap. For the (100) surface 
this gap reduction is around 0.6 eV [sij, while for the 
(110) surface we obtain a much stronger reduction, the 
surface states lie 3.3 eV below the conduction-band min- 
imum, i.e., close to the middle of the bulk GGA band 
gap, which amounts to 4.97 eV. The presence of the sur- 
face state strongly affects the screening properties with 
the consequence that the U parameter is considerably 
reduced at the (110) surface, while it remains slightly 
above the bulk value in the case of the (100) surface. 
Analogously, localized states in the gap of an insulator 
stemming from defects such as vacancies, interstitial or 
substitutional impurities at surfaces lead also to a reduc- 
tion of the surface U, but it is expected to be local, i.c, 
only the Hubbard U parameters of the defect atom and 
nearby atoms arc affected. 

Finally we discuss the frequency dependence of the sur- 
face U{uj) considering the non-magnetic state of bcc Cr 
and fee Ni. The obtained results are presented in Fig. [3] 
and compared with the bulk calculations. We see that. 



apart from the substantial reduction of the static U at 
the Cr surface, its frequency dependence is quite differ- 
ent from that of the bulk: Us{uj) increases monotonically 
with frequency, and at about 2 eV it crosses the C/b('^) 
curve towards larger values, while Ub{i^) stays almost 
constant between and 4 eV. For frequencies up to the 
plasmon frequency, Us{i^) is mostly larger than J7b(w). 
This behavior is not surprising because in the TMs the 
3d states close to the Fermi level are much more affected 
by the presence of the surface than the states at higher 
energies. Hence, with increasing frequency the polariza- 
tion function becomes less sensitive to subtleties of the 
surface electronic structure, and at high frequencies the 
surface U (ui) tends to be larger than the bulk value due 
to the dimensionality effect. 

In conclusion, by employing first-principles calcula- 
tions in conjunction with the constrained random-phase 
approximation we have determined the strength of the 
Hubbard U parameter at metal and insulator surfaces. 
We showed that U does not always increase at the sur- 
face as commonly expected. In fact, it decreases at most 
of the TM surfaces as well as insulator surfaces with pro- 
nounced surface states. We found that surface states and 
the effective band narrowing play an important role for 
the strength of the Hubbard U at surfaces. The fact 
that the U parameter can be made to increase as well as 
decrease at surfaces offers new possibilities in designing 
materials with tunable correlations. 
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